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ABSTRACT 


Pressure  drops ^  gas -liquid  discharge  ratios  at  floifing 
conditions  and  gas-liquid  test  section  ratios  were  determined  for 
air -’prater  mixtiures  flOT-ring  upifard  in  tuhes  having  inside  diameters 
of  0.630^  1.025^  1.50  and  2.50  inches  and  test  section  lengths  of 
22.88  feet.  For  constant  superficial  uater  velocities  of  0.0695_,  0,135 
0.295,  0.625,  0.869,  1.50,  2.05,  2.51,  3.47  and  7.35  feet  per  second 
the  air  rate  va^ried  from  O.OOO5II  to  O.523  cubic  feet  per  second  at 
average  test  section  condj-tions.  All  tests  vere  conducted  at  an 
average  pressure  of  3^  psia. 

At  loT-rer  water  rates  the  unit  pressure  drop  curve  was 
found  to  pa.ss  through  a  first  minimum,,  a  maximum  a.nd-  a  second  minimum. 
Using  these  inflection  points  the  curve  was  subdivided  into  four 
regimes  I,  II,  III  and  P/  characterized  by  a  decreasing,  an  increasing, 
a  decreasing,  0.nd  finally,  an  increasing  unit  pressu.re  drop,  with 
incraasing  air-water  volume  ratio.  At  higher  water  velocities  the 
maximum  and  second  minDmum  merged  a.nd  the  unit  pressure  drop  curve 
become  para.bolic  in  shape. 

The  f3.ow  patterns  were  described  using  the  pressLire  d.rop 
regimes  as  a  basis  for  their  location.  At  low  water  velocities  bubble, 
slug,  froth,  ripple,  film  and  mist  flow  are  all  present.  However, 
at  the  higher  water  velocities  froth  flow  was  not  observed. 

Tubing  diameter  had  no  effect  on  the  locus  of  the  first 
minimum  but  did  effect  the  loci  of  the  maximum  and  second  minimum. 
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The  former  \re.s  found  to  occur  at  lover  air~vater  voluiiae  ratios  as 
diameter  incres.sed  vliile  the  second  minimum  occurred  at  higher  s.ir- 
vater  volume  ratios  vith  increa^sing  diameter. 

The  pressure  drops  vere  separated  into  tvo  components^ 
the  hydrostatic  head  and  irreversibility  components^  hy  a  method 
first  suggested  by  Radford.  Superficial  friction  factors;,  ba.sed  on 
the  f loving  liquid  phase ^  and^Ldefined  by  a  Fanning  type  expression^ 
vere  calculated  and  compared  to  reduced  forms  of'..,the  superficial 
Rei^niold’s  n^ombers.  These  friction  factors  vere  fotind  to  be  as  much 
as  50^000  times  greater  than  the  friction  fa.ctors  for  single  phase 
flov.  Diameter  has  a  large  effect  on  these  friction  factors.  As 
diameter  is  increased  the  friction  factor  increases  markedly  for  the 
same  superficial  Reynold's  number. 

The  3.iquid  holdup  vas  found  to  be  dependent  upon  diameter 
and  also  upon  the  flov  pattern^  and  van  described  quantitatively  by 
a  holdup  ra.tio.  This  is  simply  a  ra.tio  betveen  the  ahr-vater  dis¬ 
charge  ratio  and  the  test  section  ratio.  The  effect  of  diameter  on 
the  holdup  ratio  vas  found  to  be  very  complex^  vith  a.  separate  curve 
existing  for  each  vater  ra.te  and  tube  diameter. 
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CliAPTM^  I 

INTRODUCTION 


The  phenomena  of  the  flow  of  gas-liquid  mixtures  ha^received 
increasing  attention  during  the  last  thirty  years.  Outstanding  examples 
of  the  flow  of  two  phase  systems  include  the  flow  of  gas-oil  mixtures 
in  the  production  tuhing  of  oil  wells  and  the  flow  of  gas -liquid  mix¬ 
tures  in  long  tube  vertical  evaporators  and  in  boiling  water  reactors. 
Design  of  these  tubing  installations  is  often  carried  out  by  trial  and 
error  methods.  A  more  efficient  sizing  of  the  tubing  requires  a  knowledge 
of  the  pressure  drop  encountered  in  the  installation.  Accurate  prediction 
of  this  pressure  drop  can  only  be  made  if  one  has  an  understanding  of 
the  mechanisms  and  principles  governing  two  phase  flow. 

Iviuch  work  along  these  lines  has  been  reported  in  the  liter¬ 
ature  but  no  complete  solution  to  the  problems  attending  vertical  two 
phase  flow  has  yet  been  presented.  This  is  undoubtedly  due  to  the 
complexity  of  this  type  of  flow  and  the  many  differences  between  single 
phase  and  two  phase  flow.  Govier^  Radford  and  Dunn  (6)  summarize  these 
differences  as  follows: 

(1)  The  flow  pattern  cannot  be  described  simply  (as  for 
single  phase  systems)  as  laminar^  transitional  or  turbulent.  In  the 
complex  flow  patterns  of  two  phase  systems  the  relative  amounts  of  the 
phases;  their  dispersion  one  in  the  other;  and  their  individLial  motions 

ti 

are  all  important. 

(2)  The  two  phase  mixture  actually  in  motion  in  a  vertical 
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tube  is  never  of  the  same  overall  composition  as  that  admitted  to  or 
released  from  the  tube^  even  under  conditions  of  steady  state.  This 
difference  in  relative  quantities  of  the  two  phases  has  been  referred 
to  under  the  headings  of  ‘slip’  and  ’holdup'.  In  general  the  fluid  of 
lower  density  tends  to  slip  pa^st  that  of  the  higher  density  so  that^ 
in  the  case  of  a  gas-liquid  system^  the  gas-liquid  ratio  in  the  flow 
tube  is  reduced  over  that  of  the  entering  or  leaving  mixtures". 

The  complexity  of  two  phase  flow^  as  compared  to  single  phase 
flow^  may  be  further  illustrated  by  considering  the  variables  involved. 
These  include  the  following: 

LIQUID  H^QPERTIES 

Density^  Pr  ^  specific  volume^  vj^. 

Viscosity,,  |j,^ 

Surface  tension,,  X 
GAS  PROPERTIES 

Density,,  or  specific  volume,,  Vq 
Viscosity,  jj,Q 
GEOliSTRICAL  PROPERTIES 

Tube  diameter,  D 
Tub  e  r ougtoe  s  s ,  o/d 
Tube  height,  X 

FLOW  RATES 

ilass  flow  rate  of  the  liquid  phase,  L 
iiass  flow  rate  of  the  gas  phase,  G 
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Previous  investigations^  in  particular  those  of  Gosline  (5) 
and  Poettraann  and  Carpenter  (ll);,  suggest  that  the  effect  of  surface 
tension  is  small.  Experimental  vork  carried  out  by  Dunn  would 
seem  to  support  this  conclusion.  Poettmann  and  Carpenter  have  sho’tm 
that  in  regions  of  turbulent  motion  gas  phase  viscosity  is  not  im¬ 
portant.  \^ile  analogy  to  single  phase  flow  necessitates  the  inclusion 
of  relative  roughness  as  a  variable;  its  effect  is  likely  submerged  by 
the  roughness'  associated  with  the  various  flow  patterns. 

V/ork  on  the  vertical  flow  of  two  ph8.se  mixtures  at  the 
University  of  Alberta  was  begun  in  19^8  by  Radford  (l2).  He  studied 
the  vertical  flow  of  air-water  mixtures  in  order  to  make  visual  ob¬ 
servations  of  the  flow  patterns  and  to  suggest  a  preliminary  method  of 
interpreting  experimental  pressure  drop  data.  This  investigation  wa.s 
continued  in  1952  by  Durm  (4).  He  further  extended  the  visual  observ¬ 
ations  taJien  by  Radford  and  used  the  method  suggested  by  Radford  in 
correlating  experimental  pressure  drop  data.  The  data  obtained  by 
Radford  and  Dunn  were  extended  and  recalculated  and  published  by  Govier; 
Radford  and  Dunn  (6). 

The  present  project  \ms  undertaken  as  a  continuation  of  the 
work  of  Dunn.  The  equipment  was  enlarged  and  improved  in  order  that 
the  effect  of  tube  diameter  on  flow  pattern;  holdup  and  pressure  drop 
which  occurred  during  vertical  two  phase  flow  could  be  investigated. 

The  equation  developed  in  the  earlier  work  is  used  to  compute  irreversi¬ 
bility  losses.  These  in  turn  are  used  to  compute  superficial  friction 
factors  (based  on  the  flowing  liquid  phase  alone)  which  are  defined 
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by  means  of  a  Fanning  type  expression.  These  friction  factors  are 
compared  to  a  reduced  form  of  the  Reynold’s  numbers ^  also  based  on  the 
flowing  liquid  phase  alone.  The  flow  patterns  are  described  using  the 
minima  and  maximum  of  the  pressure  drop  curves  as  a  basis.  These  in¬ 
flection  pointsman!  hence  the  flow  patterns^  were  found  to  be  dependent 
upon  the  superficial  vrater  velocity^  pipe  diameter  and  the  air-water 
volume  ratio^  measured  at  average  test  section  conditions. 
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CI-IAPTER  II 

REVIE  V7  OF  LITERATURE 

The  general  study  of  vertico.1  t\io  phase  flov  may  he  broken 
doim  into  an  investigation  of  the  effect  of  the  primary  variables  on 
pressure  drop^  holdup  and  flov  pattern.  One  vould  expect  that  the 
pressure  drop  could  be  correla.ted  in  terms  of  the  pertinent  variables^ 
without  a  separate  consideration  of  the  holdup  and  flow  patterns.  PIow- 
ever^  the  pressure  drop  curve  would  be  expected  to  reflect  the  change 
from  one  flow  pattern  to  another  and  the  holdvip  or  slip  velocity  to 
be  related  in  some  wa,y  to  the  flow  pattern. 

A.  PRESSURE  DROP 

Perhaps  the  most  time  and  effort  has  been  spent  in  an  attempt 
to  obtain  a  correlation  between  the  variables  involved  in  the  vertical 
flow  of  t\io  phase  mixtures  and  the  experimentally  measured  pressure 
losses.  Many  investigations  have  concluded  that  a  solution  to  this 
problem  exists  in  a  form  similar  to  that  used  for  single  phase  flow. 
This  concept  was  first  introduced  by  Uren  et  al  (l6)  and  Ilowels  (lO). 
Uren  and  his  co-workers  measured  pressure  losses  for  air-oil  mixtures 
flowing  in  a  vertical  2- inch  diameter  pipe^  4l  feet  hi^.  They  sug¬ 
gested  that  all  of  the  variables  could,  be  included  in  a  Farming  type 
expression: 

AP  .2.  fL  . (l) 

Sc  » 

where  AP  ,=  pressure  drop^  lbs. per  sq.  ft. 
f  =  friction  factor 
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L  =  length  of  section;,  feet 

p  =  density  of  the  rniicture^  Ihs.  per  cu.  ft. 

V  =  velocity  of  the  mixture^  ft.  per  sec. 

D  =  diameter  of  pipe^  ft. 

They  assumed  that  the  viscosity  of  the  rnix^xire  throughout 
the  entire  range  of  flow  conditions  coLild  he  ca.lculated  using  Poiseullle's 
equation  for  laminar  motion.  This  assumption  resulted  in  erroneous 
values  for  the  viscosity  of  the  mmctui-e.  In  presenting  their  results 
the  authors  do  not  eicplain  how  the  velocity  term  in  the,  Reynold  ’  s  number 
WS.S  calculated.  The  notion  of  slip  velocity  was  a.pparently  not  perceived 
hy  them. 

Rowels  (lO)  studied  gas  lift  data  for  flow  through  tu.hlng 
varying  from  2-l/4  to  6-l/k  inches  in  diameter.  lie  also  assumed  that 
the  Fanning  equation  could  he  used  in  calcxilating  the  pressure  drop. 
Average  values  of  the  density^  velocity  and  viscosity  of  the  mixture 
were  calculated  using  as  a  basis  the  mean  pressuire  in  the  tubing. 

The  ma-jor  difficulty  encountered  in  the  approach  of  Uren  et 
al  and  Rowels  was  in  the  determination  of  the  actual  values  of  the 
terms  density^  viscosity  and  velocity  to  be  used  in  the  calculations. 

In  fact^  no  single  density;,  velocitj^  or  viscosity  can  be  assigned  to  a 
two  phase  system. 

Martinelli  et  al  (8)  have  shown  that  four  types  of  flow 
mechanisms  can  exist  during  the  flow  of  two  phase  systems.  These  are: 

(1)  Flow  of  both  the  liquid  and  the  gas  may  be  turbulent 

(2)  Flow  of  the  liquid  may  be  turbulent  and  flow  of  the  gas 
may  be  viscouS;, 

(3)  Flow  of  the  liquid  may  be  viscous  and  flow  of  the  gas 
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may  be  tnrbulent. 

(4)  Flow  of  both  the  liquid  and  the  gas  may  be  viscous. 
Lockhart  and  Kartinelli  (7)  have  correlated  the  pressure  drop  result¬ 
ing  from  these  flow  mechanisms  in  terms  of  two  parameters^  0  and  X. 

0  is  the  ratio  of  the  pressure  drop  per  unit  of  length  during  two 
phase  flow  to  the  pressure  drop  per  unit  of  length  for  either  phase 
flowing  alone.  X  is  the  ratio  of  the  pressure  drop  per  unit  of  length 
for  liquid  flowing  alone  to  the  pressure  drop  per  unit  of  length  for 
the  gas  flowing  alone.  The  parameters  0  8.nd  X  were  evaluated  by  the 
authors  for  each  of  the  four  flow  mechanisms.  The  correlation  presented 
by  LocMiart  and  iJa^rtinelli  applies  only  to  horizontal  flow  but  suggests 
a  method  of  interpretation  that  could  be  used  for  the  irreversible 
portion  of  the  pressure  drop  in  the  vertical  flow  of  two  phase  mixtures. 

Until  recently  few  attempts  have  been  m8.de  to  separate  the 
total  pressure  drop  into  its  two  components  -  the  irreversibility  com¬ 
ponent  and  the  hydrostatic  head  component.  In  most  publications  it  has 
been  assumed  that  both  could  be  included  in  the  Fanning  expression. 

A  thermodynamic  approach  to  this  iDroblem  was  reported  by 
Relaford  (l2)  in  1948.  He  applied  a  mechanical  energy  bsAance  to  the  two 
flowing  fluids  and  developed  a  general  equation  which  permitted  the 
separation  of  the  total  pressure  drop  into  its  two  components.  This  same 
equation  was  also  used  by  Dunn  (4)  and  is  employed  here.  It  may  be 
simplified  into  two  forms^  as  shovrn  by  Govier;,  Radford  and  Dunn  (6). 

The  se  are : 
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P 


specific  voliune  of  liquid^  ft.  per  lb. 

3 

specific  volume  of  gas^  ft.  per  lb. 
pressure^  lb.  per  sq.  ft.  absoliite 
height^  feet_, 

gas-liquid  mass  ratio_,  G/L 
gas-liquid  volimie  ratio^  G  Vq 

“TTT" 


(G  =  rate  of  supply  of  liqiiid^  lbs.  per  hr.) 

(L  =  rate  of  supply  of  gas^  lbs.  per  hr.  ) 

1 

A  two  phase  friction  factor;,  fj^  j  was  defined  through  the  relation: 

/  AF\  =  2  f T.  vf.  . (it) 

Wi  Sc® 

where  g^  =  dimensional  conversion  'J^^^ctor  ft/lb.^  sec. 

=  superficial  liquid  velocity  based  on  the  total  tube 
cross  section^  ft.  per  sec. 

D  ~  tube  diameter ;  feet 

Similarly  a  two  phase  friction  factor^  f^^  was  defined 
through  the  relation: 

1  P 

. (5) 


where 


superficial  gas  velocity  based  on  the  tota.l  tube 


cross  section^  ft.  per  sec. 
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It  is  worth  noting  that  the  units  of  equation  (2)  are  feet 
of  liquid  per  foot  of  tube  and  the  units  of  equation  (3)  are  feet  of 
gas  per  foot  of  tube.  Also  the  term  Vp  / '  is  equivalent  to  the 
"submergence  ratio"  referred  to  in  articles  on  air  lift  pumps. 

Radford  (l2)  investigated  the  vertical  flow  of  air-water 
mixtiu-es  at  varying  liquid  flow  rates  in  a  1.025  inch  diameter  plastic 
pipe^  22.88  feet  in  length.  His  tests  were  carried  out  at  average 
tubing  pressures  between  16  and  23  psia.  He  separated  the  total 
pressure  drop  into  its  two  components  and  suggested  that  the  irrever¬ 
sibility  component  might  best  be  correlated  in  terms  of  a  modified 
Reynold’s  number  based  on  the  air  rs.te  alone.  Dunn  (4)  continued 
Radford’s  work  and  studied  the  vertical  flow  of  air -water  mixtures  in 
two  vertical  plastic  pipes  O.63  and  1.025  inches  in  diameter  and  22.88 
feet  in  height.  Tlie  average  tubing  pressures  studieoi  were  18  psia 
and  36  psia  at  constant  water  velocities  of  0.0695^  0.135^  0.295  and 
0.66  feet  per  second.  The  measured  pressure  drops  were  separa^ted  into 
their  two  components.  Superficial  friction  factors^  based  on  the 
flowing  gas  phase ^  and  defined  by  a  Fanning  expression^  were  plotted 
against  superficial  Reynold’s  numbers^  also  based  on  the  flowing  gas 
phase  occupying  the  entire  cross  section  of  the  tube.  No  attempt  was 
ma,de  to  determine  the  effects  of  diameter  and  gas  phase  density  on  the 
reported  friction  factors. 

Poettmann  and  Carpenter  (ll)  have  applied  a  mechanical  energy 
balance  to  the  flowing  system^,  but  in  a  much  more  approximate  manner. 
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Tliey  presented  data  comprising  a  wide  range  of  operating  variables 
obtained  from  49  flowing  and  gas  lift  wells.  Tlie  tubing  sizes  varied 
from  2  to  2-1/2  inches  in  diameter  and  ranged  in  depth  from  1000  to 
10 j 000  feet.  Tlie  tubing  pressures  varied  from  85  to  2^440  psia.  In 
attempting  to  correlate  the  energy  losses  with  the  variables  involved 
Poettmann  and  Carpenter  treated  the  flowing  fluid  mi^cture  as  a  single 
homogeneous  phase.  The  total  flowing  density  was  used  in  preference  to 
the  "in  situ"  density.  The  authors  assumed  that  the  flviid  flowing  in  the 
vertical  string  was  in  a  hi^  degree  of  tui’bulence  and  concluded  tha^t 
the  energy  loss  from  viscous  shear  was  negligible;,  and  that  the  flow 
was  independent  of  viscosity  effects.  Tlie  general  correlation  presented 
is  analagous  to  that  for  single  phase  flow  and  consists  of  a  plot  of  a 
friction  factor^  "f’4  versus  a  pseudo  Reynold’s  number.  Tlie  friction 
factor  was  defined  using  a.  Fanning  type  expression. 

Calvert  and  Williams  (2)  made  an  extensive  study  of  the 
pressure  losses  encountered  in  the  film  flow  region.  They  used  da,ta 
published  by  Radford  (l2)  and  also  conducted  further  tests  in  a  one- 
inch  diameter  copper  tube  using  water  velocities  varying  between  O.O817 
and  0.45  feet  per  second.  A  drag  coefficient  was  used  to  compute  the 
pressure  losses  due  to  profile  drag.  This  was  defined  by  the  follow¬ 
ing  equation^  derived  by  the  authors: 


A?. 


%  Vq  Pq  Ap 


(6) 


lihere 


pressure  drop  due  to  drag^  lb.  per  sq.  ft. 


drag  coefficien- 
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Vq  =  superficial  velocity  of  gas  phase  hased  on  total 
cross  sectional  area^  ft.  per  sec. 

=  gas  density ;  Ih.  per  cu.  ft. 

Ap  -  projected  area  for  drag;  sq.  ft.  per  sq.  ft. 
inside  pipe  area. 

L  =  length  of  tube;  feet 

D  =  diameter  of  tube;  feet 

g  =  [dimensional  conversion  factor  lb.  ft.  per  Ib.^  sec. 

The  value  for  vas  determined  experimentally  for  the  conditions 
studied  and  presented  as  a  function  of  film  thickness;  the  mass  flov 
rate  ofuthe  gas  phase  and  the  projected  area  for  drag.  The  authors 
postulated  a  method  of  determining  the  total  pressure  drop  by  adding 
the  value  of  the  skin  friction  losses  to  the  profile  drag  losses.  The 
la.tter  could  be  evaluated  using  equation  (6).  In  general  this  method 
of  predicting  pressure  drop  is  not  applicab3.e  to  the  entire  range  of  air- 
uater  ratios  beca.use  of  the  V8.rying  and  complex  nature  of  the  flou  patterns 
involved . 


B.  FLai  PATTERIIS 

Basically  the  flov  of  gas-liquid  mixtures  may  be  said  to  be 
one  of  four  type: 

(1)  A  dispersed  phase  of  gas  bubbles  throughout  the  liquid 

(2)  A  dispersed  phase  of  liquid  droplets  throughout  the  gas. 

(3)  Alternate  slugs  of  gas  and  liquid 

(4)  A  core  of  gas  flowing  through  an  annular  ring  of  liquid 
The  first  investigation  of  flow  patterns  occurring  in  the 
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vertical  flow  of  two  phase  mixtures  was  published  by  Cromer  and  Hunting- 
ton  (3).  They  have  also  been  described  by  Gosline  (5)^  Bergelin  (l); 
Stovall  and  Huntingdon  (15)  Calvert  and  Williams  (2)  Radford  (l2)^ 

Dunn  (4)  Stein  et  al  (l4)^  and  Govier^  Radford  and  Dunn  (6). 

Cromer  and  Huntington  made  visual  studies  of  the  flow  of  air- 
water  mixtures  in  a  2-inch  diameter  vertical  pipe  approximately  90  feet 
in  length.  Observation  sections  were  placed  at  l4  foot  intervals  along 
the  pipe  and  the  flowing  mixture  was  photographed  at  these  points. 

Most  of  the  descriptions  of  flow  patterns  were  ba^sed  upon 
"efficiency"  or  "yield"  considerations.  Lader  experimenters^  notably 
Radford  (l2)^  Dunn  (4)^  Calvert  and  Williams  (2)  and  Govier^  Radford 
and  Dunn  (6)  have  described  the  flow  patterns  in  conjunction  with  pressure 
drop  measurements.  The  latter  subdivided  each  pressure  drop  curve  into 
regimes;,  and  on  the  basis  of  these  have  described  the  flow  patterns. 

This  method  of  description  would  seem  to  be  more  precise  than  those 
used  by  Cromer  and  Huntington  and  other  earlier  investigators^  and  is  fol¬ 
lowed  here. 

The  flow  padterns  occurring  in  the  vertical  flow  of  air-water 
mixtures  can  best  be  described  with  the  aid  of  the  unit  pressure  drop- 
air-water  volume  ratio  relationships  such  as  are  presented  later  as 
Figures  8  to  25.  These  pressure  drop  curves  may  be  divided  into  four 
regimes^  III  and  IV^  using  the  maximum  and  minim-um  points  as  the 

basis  for  the  division.  These  regimes  a.re  characterized  respectively  by 
a  decreasing^  an  increasing,  a  decreasing  and  finally  an  increasing 
pressure  drop,  with  increasing  air-water  volume  ratio.  In  attempting  to 
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describe  the  flow  patterns  one  is  confronted  with  the  problem  that  the 
regimes  and  corresponding  flow  patterns  vary  slightly  with  changing  welter 
rates.  For  this  reason  it  is  convenient  to  treat  the  high  and  low  water 
rates  separately. 

(l)  Low  VJater  Rates 

At  very  low  air-water  discharge  ra^tios  the  flow  consists  of 
bubble  dispersed  gas.  This  type  of  flow  occurs  in  a  very  narrow  region 
at  the  beginning  of  regime  and  is  illustrated  in  Fignre  1.  As  the  adr 
rate  increases  these  bubbles  become  3.arger  in  size  and.  coa.lesce  to  form 
bullet- shaped  slugs.  This  type  of  flow^  shown  in  Figure  2^  consists  of 
alternate  slugs  of  air^  su.rrounded  by  a  thin  water  annulus^  and  bubble 
dispersed  ga.s.  After  a  slug  passes  any  given  section  of  the  tube  the 
water  appears  to  fall  d.o^-Tnwards  a  short  distance  to  form  a  region  of 
bubble  d-ispersed  gas  beneath  it.  Slug  flow  continues  beyond  the  first 
minimuju  and  about  half^ray  into  regime  II.  There  are^  however^  at  least 
two  differences  in  slug  flow  between  regimes  I  and  II.  In  regime  II 

(i)  the  length  and  velocity  of  the  individual  slugs  has 
increased  and  the  slug  spacing  has  decreased^ 

(ii)  the  inner  surface  of  the  liquid  annulus  of  the  slug 
has  become  wary  in  appearance. 

About  midway  between  the  first  minimum  and  the  maximun  the 
slugs  are  replaced  by  a  frothy  pattern  which  extend.s  a^long  the  entire 
leng-th  of  the  tube.  This  type  of  flow  is  characterized,  by  the  extreme 
tui-bulence  of  the  liquid  phase ^  8.nd  is,  illiistrated  in  Figure  3*  "^^e 

region  of  froth  flow  extends  from  the  middle  of  regime  II  a.lmost  to  the 
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Figure  1.  Bubble  Flow 
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Figure  2.  Slug  Flow 
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Figure  3»  Froth  F1o\t  At  A  Loir  Liquid  Velocity 
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end  of  regime  III.  As  the  second  minimum  is  approached  a  gra.dation 
from  froth  to  ripple  flow  is  noticed.  This  latter  type  shown  in 
Figure  4  is  characterized  hy  the  ripply  or  wavy  nature  of  the  liquid 
phase . 

With  a  further  increase  in  the  air  rate  the  ajnount  of  dis¬ 
turbance  in  the  liquid  phase  decreases  steadily^  and  about  the  begin¬ 
ning  of  regime  IV  film  flow  is  formed.  Figure  5  illustrates  this  type 
of  flow  pattern.  Eventually  at  some  very  hi^  air  rate  waiter  droplets 
would  be  captured  from  the  air-liquid  interface  resulting  in  mist  flow. 

(2)  High  Water  Rates 

At  low  air  rates  the  flow  patterns  a.re  identica-l  to  those 
previously  described^  i.e.^  the  wa,ter  rate  has  no  appreciable  effect  on 
bubble  dispersed  gas  or  slug  flow. 

Once  the  region  corresponding  to  froth  flow  is  entered  the 
pattern  seems  to  change  slightly  from  thab  occurring  at  lower  vatev 
rates.  The  liquid  phase,,  as  before^  is  characterized  by  a  frothy  pattern 
which  extends  along  the  entire  length  of  the  tube.  However^  fast  mov¬ 
ing  air  slugs  travel  intermittently  through  the  tube;,  interrupting  the 
continuity  of  the  liquid  phase.  Figure  6  shows  this  flow  pattern. 

As  more  air  is  admitted  to  the  system  these  slugs  become 
faster  moving  and  more  frequent  in  occurrence.  At  some  sufficiently 
high  air  rate  ripple  flov^  will  result.  The  remaining  two  types  of  flow^ 
film  and  mist^  are  unchanged  at  high  liquid  velocities. 

To  summarize ;  the  essential  difference  in  flow  patterns  be- 
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Figure  4 


Ripple  F1o^7  At  A  ]jD\r  Liquid  Velocity 
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Figure  5*  Fi3jn  Flovr 


Figure  6.  Froth  Flow  At  A  High  Liquid  Velocity 
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tween  low  and  high  water  rates  is  the  mechanism  of  transition  from  slug 
to  film  flow.  In  the  case  of  high  water  rates  the  froth  flow  region  is 
characterized  hy  the  presence  of  air  slugs.  One  shoiild  remember  that  the 
boundaries  describing  the  transitions  from  one  flow  pattern  to  another 
are  not  too  distinct  and  their  positions  on  the  pressure  drop  curve  are 
someirhat  arbitrary. 

C.  HOLDUP 

Moore  and  Wilde  (9)  were  among  the  first  to  recognize  the 
presence  of  slippage  in  the  vertical  flow  of  two  phase  mi^rbures  and  to 
make  a  quantitative  study  of  the  variables  affecting  it.  They  conducted 
extensive  tests  of  slippa^ge  in  1^  2^  2-1/2;,  3  s-nd  4  inch  diameter  pipes 
approximately  67  feet  high.  Five  liquids  were  used:  water^  kerosene^ 
heavy  gas  oil^  light  lubricating  oil  and  heavy  lubricating  oil.  They 
believed^  that  kno^ring  the  fraction  of  liquid  in  the  pipe^  it  would  be 
possible  to  ca^lculate  the  liquid  hydrostatic  head.  However^  as  Radford 
(12)  demonstrated^  this  assumption  was  erroneous.  Because  of  this  their 
interpretive  work  is  of  little  fundamental  value^  but  their  discussion 
does  give  one  a  clearer  under sta^nding  of  holdup. 

Gosline  (5)  used  data  on  the  study  of  single  bubble  motion 
as  a  beginning  to  the  investigation  of  slippage.  From  theoretical  con¬ 
siderations  he  derived  an  equation  for  calculating  the  mean  velocity 
of  the  gas  using  experimentally  measurable  quantities.  He  a.lso  investi¬ 
gated  the  flow  of  air -liquid  mixtures  in  a  one -inch  diameter  vertical 
pipe^  using  as  liquids^  water  and  two  colorless  petroleum,  oils.  Using 
these  data^  a  drag  coefficient  (based  on  an  equation  derived  in  the  paper) 
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was  plotted  versus  the  air  core  velocity.  This  coefficient  diminished 
rapidly  with  increasing  air-velocity. 

In  their  studies  of  slip  velocity  Stein  et  al  (4)  used  a 
modification  of  an  equation  developed  hy  Gosline.  The  data  were 
obtained  in  gas -liquid  systems  which  included  combinations  of  air^ 
propane  and  natural  gas  with  water ^ .  lubricating  oils  and  crude  oils. 
The  tests  were  conducted  in  2,  3  s-nd  4  inch  diameter  tubes.  The 
average  gas  phase  density  varied  from  0.0426  to  O.7066  pounds  per 
cubic  foot.  The  authors  concluded  that^  of  all  the  variables  studied^ 
gas  rate  and  pipe  size  had  the  greatest  effect  on  the  slip  velocity. 

Radford  (l2)^  Dunn  (4)^  and  Govier^  Radford  and  Dunn  (6) 
have  reported  investigations  of  holdup  occurring  in  the  vertical  flow 
of  air-water  mixtures.  Their  stiidies  were  conducted  using  O.63O  and 
1.025  inch  diameter  pla.stic  tubes  at  average  gas  phase  pressures  of 
18  and  36  psia.  The  concept  of  holdup  ratio  was  intioduced  by  Govier^ 
Radford  and  Dunn  (6).  It  is  simply  the  ratio  between  the  flowing 
air-water  ratio  and  the  *'in  situ”  or  test  section  air-water  ratio. 

The  term  holdup  ratio  is  used  as  a  quantitative  measure 
of  the  development^  within  the  flow  section^  of  a  concentration  of 
the  denser  phase  greater  than  that  in  the  supply  mixture.  Slip 
velocity  is  defined  as  the  difference  in  the  average  velocities  of  the 
gas  and  liquid  phases^  both  measLired  at  average  flowing  conditions. 
Since  both  quantities  are  measures  of  the  same  phenomenon  it  follows 
that  they  are  interrelated.  It  may  be  easily  sho^.-m  (6)  that  this 
relationship  is  given  by  the  following  equation: 
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where 


(7) 


holdup  ratio  (defined  above) 
average  velocity  of  gas  phase 
average  velocity  of  liquid  phase  in  the  same 
units  as 


ft-  ■ 
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CHAPTER  III 

EXPERIMENTAL  EQUIPMENT 

The  equipment  is  illustrated  schematically  in  Figure  7. 

It  is  essentially  the  same  as  that  used  hy  Radford  and  Dunn.  It 
consists  of  an  air  supply  and  metering  system,,  a  water  supply  and 
metering  system,,  four  test  sections  of  diameters  O.63,,  1.025^  I.50 
and  2.50  inches,,  and  a  manometer  system  for  measuring  the  pressure 
drop. 

The  air  was  supplied  from  a  central  compressor  through  a 
reducing  valve  x/hich  maintained  a  constant  pressure  of  75  psig.  to 
the  orifice  metering  system.  The  volume  of  air  admitted  to  the  hase 
of  the  equipment,,  through  the  metering  system,,  X7as  controlled  manu¬ 
ally.  Rough  control  was  effected  hy  a  valve  in  the  one-inch  supply 
line,,  while  a  needle  valve  in  a  1/4  inch  hy-pass  line  made  finer 
control  possible.  Water  was  supplied  through  two  70  gallon  per 
minute  centrifugal  pumps,,  connected  in  series,,  having  an  available 
hea.d  of  220  psig.  The  water  was  metered  using  one  of  four  rotameters 
having  maximum  capacities  of  O.OOO272,,  0.00247,,  0.0011  and  O.O7  cubic 
feet  per  second. 

Preliminary  tests  indicated  that  the  type  of  mixing 
section  had  no  appreciable  effect  upon  tvo  phase  flow  phenomena.  For 
this  reason  the  air  and  water  x-zere  brought  together  using  a  simple 
tee.  The  four  test  sections  xrere  connected  in  parallel,,  each  section 
having  its  ox7n  mixing  tee.  The  O.63  and  1.025  inch  tubes  X'/ere  both 
of  extruded  cellulose  acetate  butyrate  plastic^  the  I.50  inch  tube 
X7as  standard  seamless  copper  water  tubing  and  the  2.5O  inch  tube  was 


¥ 


O 


.;^c 


:  O'  c'lC' 


i  >: 


,  w>i«* 


:£:;, :  o  ';:o' 


. .  I.. 


0 . 


O' 


D.’> 


D.'. 


,'0 


'/■  .  O' 'O.:;': 


(  :  .  .  - ,  .  o  . 


i'O.J 


.  .0  A.  0:0'::,  ,  t/T  c;  ...  o 

O:., ‘/.o^ ':  r ‘'y  i:<i  d  o.  O' ■  ;rj:dci‘: ' 
>".o  :  .0:;  o^:w'o£i|;r  ;;:v0/-  ■  . '  '  ;  .  ■ 

'c'  oo'J'  ,, t  o,..  :J'.  ,  .  ori^-'  ,. 

'-i'  OM  ,,  ...  okJ’  o)/,Jo;.'ov;'  3vO  r r  Ov^  ;■■;■■  CovL’’;)-"  ',  .  J'.'o,'  -’v 

'  ^  '■  ,  ^  E®  OO'. 


27 


Figure  7 


SCHEMATIC  FLOW  DIAGRAM  OF  EQUIPMENT 
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Commercial,  carbon  black  loaded,  extruded  polythene.  The  lower  end 
of  each  pipe  was  connected  to  the  mixing  tee  with  an  appropriate 
length  of  galvanized  pipe  and  a  full  flow  plug  valve.  Similarly 
the  top  of  each  tube  was  connected  to  the  discharge  piping  through 
another  full  flow  plug  valve.  These  valves  operated  with  a  90  degree 
turn  of  the  handle.  Tlie  top  and  bottom  handles  of  the  two  valves  for 
each  tube  were  connected  with  a  1  inch  aluminum  angle  iron  enabling 
the  operator  to  close  both  simultaneously.  Thus  the  flowing  mixture 
could  be  trapped  in  the  test  section  and  readings  talien  of  the  "in 
situ"  ratio  of  air  to  water. 

Three  pressure  taps  were  made  in  each  tube.  The  bottom 
tap  was  pierced  approximately  8  feet  above  the  plug  valve  and  upper¬ 
most  tap  22.88  feet  above  it.  The  third  pressure  'tap  was  placed 
midway  between  the  other  two.  In  the  case  of  the  two  bu:t^''rate 
tubes  the  taps  were  made  from  blocks  of  lucite  circularly  machined 
to  fit  the  tube,  and  sealed  with  chloroform.  Tliese  taps  were  then 
bored,  carefully  reamed  and  counter  bored,  and  tapped  for  standard 
1/8  inch  copper  fittings.  The  taps  for  the  I.50  inch  tube  were  made 
by  drilling  holes  for  standard  1/8  incn  copper  fittings  which  were 
then  silver  soldered  in  place.  To  tap  the  2.5O  inch  tube  a  3/16  inch 
hole  was  drilled  through  the  tube  and  carefully  reamed.  A  l/h  inch 
pol^dihene  tube  was  then  heat  welded  over  the  hole  on  to  the  pipe. 
Connected  to  the  discharge  end  of  the  test  sections  wqs  a  large 
air-water  separator  which  permitted  water  to  be  discharged  to  waste 
and  air  to  be  discharged  tlirou^i  a  pneumatic  control  valve  to  the 
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atmosphere . 

Tiv’-o  30  inch  Meriam  manometers^  one  filled  with  water  and 
the  other  with  mercury^  were  used  to  measure  the  pressure  drop  across 
the  air  orifice  plates.  The  pressure  upstream  of  the  orifice  plate 
was  measured  using  a  0  to  100  psi.  gauge.  Connected  in  parallel  to 
the  mid  point  pressure  tap  were  a  60  inch  Meriam  manometer  and  a 
0  to  25  psi.  gauge^  Brown  pressure  recorder  a-nd  controller.  This 
controller  actuated  the  pneuma/bic  valve  located  in  the  air  discharge 
line;  and  effected  su.tomatic  control  of  the  mid  point  pressure.  Tv/o 
60  inch  I'leriam  manometers;  one  filled  with  mercury  and  the  other 
filled  with  water;  were  used  to  measure  the  pressure  differential 
Between  the  upper  and  lower  pressure  taps.  These  manometers  were 
connected  in  parallel  tlu-ough  plastic  separators  connected  directly 
to  the  pressure  taps.  Thermometers  were  provided  to  measure  the  in¬ 
let  temperatures  of  the  air  and  water. 
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CI-IAPTER  IV 

EXPERIMENTAL  PROCEDURE 

All  of  the  tests  reported  here  were  conducted  at  a  mid -point 
pressure  of  36  psia^  maintained  to  within  0.2  psia.  Tests  were  carried 
out  at  five  or  more  constant  water  rates  in  each  tube;  with  essentially 
the  S8jne  procedure  being  used  for  each  test. 

A  test  was  begun  by  opening  the  plug  valves  and  closing  all 
of  the  pressure  taps.  Water  was  admitted  to  the  test  section;  at  a  rate 
which  was  slightly  hi^er  than  that  desired;  th-rough  the  appropriate 
rotameter.  Air  was  allowed  to  enter  the  system  at  the  desired  rate 
througli  the  necessary  size  of  orifice  plate  and  the  mid-point  pressure 
tap  was  opened.  V/hen  the  mid-point  pressure  had  stabilized  at  36 
psia  the  water  rate  was  readjusted  to  the  correct  value  and  the  flow 
of  water;  discharging  to  the  sewer  from  the  large  air-water  separator; 
was  set.  Once  steady  state  conditions  vrere  reached  the  upper  and  lower 
pressure  ta.ps  were  oi^ened  and  readings  talien  of  the  water  rate;  air 
rate  and  the  pressure  drop.  IRien  these  had  been  recorded  the  upper  and 
lower  pressure  taps  were  closed  and  both  plug  valves  closed  simultaneous¬ 
ly;  enabling  the  operator  to  trap  the  flowing  mixture  between  them. 

After  a  short  period  of  time  had  been  allowed  for  the  water  to  drain 
down  the  sides  of  the  tube;  the  relative  volumes  of  air  end  water  were 
determined.  The  entire  procedure  was  then  repeated  at  other  air  rates. 

Ceneful  visual  studies  were  carried  out  in  both  the  3./2  inch 
and  the  1  inch  tube  and  several  series  of  photographs  were  taken.  These 
photographs  have  been  published  by  Govier;  Radford  and  Dunn.  Radford 
and  Dunn  also  presented  similar  photographs  in  their  theses. 
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CI-IAPTER  V 

EXPERIMENTAL  RESULTS 


The  hasic  two  phase  flow  data  are  presented  in  Tables  AI 

and  All  of  Appendix  A,  Tables  Cl  and  CII  of  Appendix  C  give  the  data 

for  the  calibrations  of  the  water  rotameters  and  of  the  air  orifice 

plates.  The  results^  calculated  from  the  data  given  in  Appendix  A, 

are  reported  in  Appendix  B.  The  friction  factors^  f^^  are  included 

L 

in  this  Appendix. 

Tlie  basic  data  given  in  Appendix  A  for  each  of  the  constant 
water  rates  in  each  tube  are  summarized  in  Figures  8  to  25.  These 
show  the  relationships  between  the  total  unit  pressure  drop  and  the 
air -water  volume  ratio  (at  flowing  conditions)  and  betv/een  the  holdup 
ratio  and  the  air -water  volume  ratio.  There  is  a  marked  similarity  in 
the  general  shape  of  the  pressure  drop  and  holdup  curves  to  those 
presented  by  Radford,,  Bunn  a,nd  Govier^  Radford  and  Dunn. 

The  bounde^ries  of  the  pressure  drop  regimes  are  shoim  on 
each  curve.  The  method  used  for  locating  each  regime  is  the  same  as 
outlined  by  Govier,,  Radford  and  Dunn.  These  boundaries  are  used  as 
a  basis  of  describing  the  flow  patterns  and  also  in  defining  the 
effects  of  diameter  on  flow  pattern. 

At  lower  water  velocities  the  existence  of  the  minima  and 
maxi murn  points  in  the  pressure  drop  curve  for  all  diameters  is  appar¬ 
ent.  At  hitler  water  velocities  the  ma:cimum  merges  with  the  second 
minimum  and^  in  fact^  both  disappear  at  a  sufficiently  large  super- 
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ficial  water  velocity.  This  merging  effect  does  not  occur  at  the  same 
liquid  rate  in  each  diameter.  In  the  O.63O  inch  diameter  tube  the 
mascimum  tends  to  blend  into  the  second  minimum  at  a  superficial  water 
velocity  of  0.86  feet  per  second.  The  corresponding  mergence  occurs 
at  a  superficial  water  velocity  of  2.5  feet  per  second  (Figure  I3)  in 
the  1.025  inch  diameter  tube.  It  v/'as  not  possible  to  determine  experi¬ 
mentally  the  point  of  disappearance  of  the  second  minimum  in  the  I.50 
inch  and  2.5O  inch  diameter  tubes. 

As  is  the  case  for  the  pressure  drop  curves,  the  general 
diape  of  the  holdup  curves  is  the  same  for  all  diameters.  A  point  of  in¬ 
flection  occurs  about  the  middle  of  regime  I.  This  inflection  point 
occurs  at  all  water  velocities  a.nd  in  all  tube  diameters,  although  its 
presence  is  most  easily  detected  at  low  liquid  rates.  Beyond  this  in¬ 
flection  point  the  holdup  curves  for  the  range  of  conditions  studied 
may  eidiibit  one  of  three  characteristics: 

(1)  The  holdiip  ratio  may  increase  approximately  linearly 
to  a  maximum  value  and  then  decrease, 

(2)  The  holdup  ratio  may  increase  approximately  linearly 
to  a  maximum  value  and  then  remain  nearly  constant, 

(3)  The  holdup  ratio  may  increase  uniformly  exhibiting 
no  abrupt  change  of  slope. 

It  should  be  remembered  that  when  either  pure  air  or  pure 
water  is  the  flowing  system  the  value  of  the  holdup  ratio  will  be  unity 
as  shown  by  equation  7.  This  suggests  that  a  holdup  curve  illustrating 
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a  ran^L-e  of  air-water  ratios  from  zero  to  infinity  must  e^chibit  a 
maximi.im  at  some  point. 

The  value  of  the  holdup  ratio  increases  as  tube  diameter 
increases  for  low  water  velocities  and  constant  air-water  volume 
ratios.  This^  however^  does  not  apply  to  the  highest  water  veloc¬ 
ities.  Here  the  holdup  ratio  for  a  larger  diameter  tube  may  be 
equal  to  or  less  than  the  holdup  ratio  for  a  smaller  diameter  tube 
at  the  same  air -water  volume  ratio.  This  complexity  is  further 
illustrated  by  Figures  3^  anl  37  described  later. 
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CPIAFPER  VI 

INTERPRETATION  OF  RESULTS 
A.  EFFECT  OF  DIAIV'IETER  ON  IRESSURE  DROP 
Tlie  total  unit  pressui-e  drop  was  separated  into  its  two 
components  -  the  hydrostatic  head  and  irreversibility  component^  by 
utilizing  equation  (2).  Using  equation  (4)  and  the  value  of  the 
irreversibility  component;  the  superficial  friction  factor  was  deter¬ 
mined.  The  numerical  values  of  this  friction  factor  are  given  in 
Appendix  B.  Figure  26  shows  the  relationship  of  this  friction  factor ; 
f  ,  to  the  air-water  volume  ratiO;  at  a  constant  velocity  for  each  of 

Xj 

the  four  tube  diameters.  As  can  be  seen  from  this  figure  the  friction 
factor  increases  with  increasing  diameter  for  typical  medium  values  of 
the  superficial  water  velocity  and  air-water  volume  ratio.  Figure  26 
clear Ij^  illustrates  the  effect  of  diameter  on  the  superficial  friction 
factor . 

Figure  27  was  prepared  from  a  complete  set  of  curves  similar 

T 

to  those  of  Figure  26.  The  former  is  a  logarithmic  plot  of  f  versus 

Li 

3*5  3 .  ^ 

DVj^;  showing  lines  of  constant  VqD  .  The  parameter  was 

determined  by  trial  and  serves  to  correlate  the  data  for  s.ll  diameters 
and  air  velocities.  It  should  be  noted  that  the  quantity  DV  iS;  in 
effect;  a  reduced,  form  of  the  superficial  Reynold's  number;  DVj^  ; 

Toyir 

unchanged  in  all  tests. 

The  conventional  single  phase  friction  factor  line  is  also 
shown  on  Figure  27  for  reference  purposes.  The  lines  representing  the 
two  phase  condition  ane  seen  to  be  well  above  and  to  the  right  of  the 
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single  phase  line^  particularly  at  low  Reynold’s  numhers.  At  low  values 
of  and  higher  values  of  DVj^  the  two  phase  lines  approach  the  single 

phase  line  and  hecome  more  nearly  parallel  to  it.  This  is  to  he  expected 
since  when  the  flowing  fluid  is  water,  i.  e.  Vq  =  0,  the  tv7o  phase 

friction  factor  must  hy  definition  he  equal  to  the  single  phase  friction 
factor . 

EFFECT  OF  DIAMETER  ON  FLOW  PATTERN 
As  pointed  out  previously  the  flow  patterns  may  he  described  in 
terms  of  the  pressure  drop  regimes.  These  regimes  are  defined  hy  the 
first  and  second  minima  and  the  first  maximum  of  the  pressure  drop  curves. 
Hence  the  effect  of  diameter  on  flow  pattern  can  hest  he  described  hy 
considering  its  effect  on  the  loci  of  the  maximum  and  minima.  The  effect 
of  diameter  will  he  reflected,  both  in  the  ordinate  (imit  pressure  drop) 
and  the  abscissa  (air-water  volume  ratio)  for  each  of  these  loci. 

Figures  28  and  29  show  the  loci  of  the  first  minimum  and 
maximum  respectively,  as  functions  of  unit  pressure  drop  and  super¬ 
ficial  water  velocity  for  lines  of  constant  diameter. 

It  is  significant  that: 

(1)  At  lower  superficial  water  velocities  the  pressure  drop 
occurring  in  a  larger  diameter  tube  is  gxeater  than  that  occurring 

in  a  smaller  diameter  tube. 

(2)  At  the  hi^er  superficial  water  velocities  the  reverse  is 
true,  i.  e.  as  tube  diameter  decreases  the  pressure  drop  increases. 

Figure  30  shows  the  locus  of  the  second  minimum  of  the 
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pressure  drop  curve  as  a  function  of  unit  pressure  drop  and  superficial 
water  velocity  for  lines  of  constant  diameter.  However^  it  will  be 
seen  that  the  pressure  drop  at  constant  superficial  water  velocity  in¬ 
creases  with  increasing  diameter.  Because  of  this  relative  simplicity 
the  curves  of  Figure  30  uere  reduced  to  a  single  relationship^  as  shown 
by  Figure  3I.  The  latter  shows  the  locus  of  the  second  minimum  as  a 

function  of  unit  pressure  drop  and  the  ratio  The  quantity 

I  0  8s 

V^/D  '  '  was  determined  by  tria,l  a.nd  serves  to  correlate  the  data 

for  all  water  velocities  and  diameters.  Tliis  rela.tionship  is  a  straight 

line  with  a  positive  slope  of  1.7. 

I.Tiile  the  above  figures  are  of  some  interest  they  do  not  lend 
themselves  to  an  easy  description  of  the  flow  patterns.  From  this  pdint 
of  view  Figures  32^  33;  3^  and  35  are  certainly  more  useful.  Figure  32 
shows  the  loci  of  the  first  and  second  minima  and  the  maximum  of  the 
pressure  drop  curve  as  a  function  of  superficial  water  velocity  and 
air -waiter  volume  ratio.  The  locus  of  the  first  pressure  drop  minimum^ 
i.e.;  the  boundary  between  regimes  I  and  II  is  a  strai^it  line  on 
logaritlirnic  paper  of  slope  equal  to  -1.2.  The  locus  of  the  maximum  for 
each  diameter  is  rou^ly  parallel  to  this  but  occurs  at  higher  air-water 
volume  ratios.  Each  diameter  is  represented  by  a  separate  curve  with 
the  curves  for  larger  diameters  occurring  at  lower  air -water  volume 
ra-tios  than  those  for  the  smaller  diameters.  The  loci  of  the  minima 
a.re  roughly  parallel  to  that  of  the  first  minimum  end  maximum  points^ 
but  are  displaced  towards  hi^ier  air-water  volume  ratios.  Again  each 
diameter  is  represented  by  a  separate  cur^/e  with  the  curves  for  la.rger 
diameters  occurring  am  higher  air-water  volume  ra.tios  than  those  for 


small  diameters. 
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Some  of  the  maximum  loci  hegin  to  approach  the  loci  of  the  second  minima 
at  large  values  of  the  superficial  water  velocity.  This  merging  effect 
"becomes  less  noticeable  in  the  larger  diameter  tubes  and  in  fact  the 
locus  of  the  m8:>cimum  for  the  2.5O  inch  diameter  tube  is  a  straight  line. 

The  data  represented  by  Figure  32  would  suggest  that  the  maximum  of  the 
pressure  drop  curve  will  not  exist  in  a  very  large  diameter  tube.  This  is 
illustrated  by  the  gradual  shifting  of  the  locus  of  the  maximum  to  that  of 
the  first  minimum. 

Figures  33^  3^  35  show  the  loci  of  the  inflection  points 

as  functions  of  superficial  waiter  velocity  and  air-water  volume  ratiO; 
for  each  of  the  diameters.  The  first  of  these  gives  the  loci  of  the  first 
minimum  as  a  function  of  superficial  water  velocity  and  air  water  volume 
ratio.  It  is  a  straight  line  on  logarithmic  paper  of  slope  equal  to  -1.2. 
Diameter  has  no  appreciable  effect  on  the  position  of  the  first  minimum. 

The  approximate  boundaries  between  bubble  and  slug  flow  and  between  slug 
and  froth  flow  are  sho^m  as  dashed  lines.  It  will  be  noted  that  the 
point  of  transition  between  slug  flow  and  froth  flow  is  dependent  to  a 
small  extent  upon  diameter.  This  transition  occurs  at  sli^tly  smaller 
air -water  volume  ratios  as  the  tube  diameter  increases^  for  any  parti¬ 
cular  value  of  the  superficial  water  velocity.  In  all  cases  the  transition 
occurs  at  an  air-water  vol-ume  ratio  greater  than  that  corresponding  to  the 
locus  of  the  first  minimum.  The  boundaries  between  slug  and  froth  flow 
are  represented^  in  the  typical  case,  by  straight  lines  roughly  parallel 
to  the  locus  of  the  first  minimum  points.  However,  the  data  for  the 
smaller  diameters  deviates  from  this  linearity  and  converges  towards  the 
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line  of  the  maximmn  points.  Tliis  convergence  corresponds  to  the  dis¬ 
appearance  of  the  second  minimum  discussed  previously. 

Figure  3^  shows  the  loci  of  the  pressure  drop  majcimum  as 

1/3 

function  of  the  superficial  water  velocity  and  the  product  R_^D  ; 
which  was  determined  hy  trial  and  serves  to  correlate  the  data  for 
all  diameters.  This  relationship  is  a  straight  line  of  negative  slope 
on  logarithmic  paper.  The  smaller  diameters  are  not  represented  hy  a 
straight  line  on  logaritlimic  paper ^  hut  the  loci  converges  towards  that 
of  the  second  minimimi.  This  convergence  first  occurs  at  superficial 
water  velocities  of  0.30^  0.60  and  O.9O  in  the  0.63^  1.025  and  I.50 
inch  diameter  tuhes^  respectively.  The  boundaries  between  slug  and 
froth  flow  and  hetireen  froth  and  ripple  flow  are  sho^ra  as  light  dashed 
lines.  The  latter  boundary  is  a  straight  line  which  converges  slightly 
towards  the  locus  of  the  pressure  dop  maximum  at  higher  superficial 
water  velocities. 

Figure  35  gives  the  loci  of  the  second  minimum  of  the  pressure 
drop  curve  as  a  function  of  the  superficial  water  velocity  and  the  ratio 

■  1/3 

R^/D  .  Tlie  latter  was  determined  by  trial  and  serves  to  correlate 
the  date  for  all  diameters.  This  locus  is  also  a  straight  line  of 
negative  slope  on  logarithmic  paper.  The  boundaries  between  froth  and 
ripple  flow  and  between  ripple  and  film  flow  are  shoim  as  light  dashed 
lines  which  are  roughly  parallel  to  that  determining  the  locus  of  the 
minimum.  The  region  of  ripple  flow  is  seen  to  be  a  narrow  band  compris¬ 
ing  the  end  of  regime  III  and  the  beginning  of  regime  IV. 

G.  EFFECT  OF  DIAIvIETER  OR  EOIDUP 

Figures  36  and  37  show  the  holdup  ratio  as  functions  of  the 
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superficial  vater  velocity  and  the  terms  ^  and  11^/d  respectively. 

V3  ,  1/3 

The  parameters  R^D  and  were  determined  hy  trial  and  serves  to 

correlate  the  data  for  all  diameters.  Figure  3^  is  a  plot  of  holdup 
ratio  versus  superficial  water  velocity  showing  lines  of  constant  R^D  ^  . 

The  locus  of  the  maximum  of  the  pressure  drop  curve  is  also  sho\7n.  It 
will  he  noted  that  this  relationship  is  a  series  of  smooth^  and  fairly 
flat  curves  in  the  region  of  slug  and  froth  flow.  Plowever^  once  the 
region  of  film  flow  is  approached  the  relationship  hrealis  dam  and  a 
separate  curve  is  obtained  for  each  diameter.  Tliis  region  is  seen  to  be 
to  the  right  of  the  locus  of  the  ma:cimum  of  the  pressure  drop  curve. 

Figure  37  shows  the  corresponding  relationship  for  regions  of 
film  flow.  It  is  a  plot  of  holdup  ratio  versus  superficial  water  velocity 

/  3 

for  lines  of  constant  R  /D  .  The  locus  of  the  second  minimimi  is  also 

V 

shorn.  It  is  apparent  that  this  relationship  is  only  valid  in  regions 

of  film  flow  ajid  possibly  ripple  flew.  Once  the  region  of  froth  flow 

is  entered  a  separate  curve  is  obtained  for  each  diameter.  In  general 

this  relationship  is  a,  series  of  straight  lines  of  fairly  steep  negative 

slope .  It  should  be  noted  that  the  holdup  ratios  at  value  of  Pv^/D  ‘  ^ 

1/3 

of  600  are  lower  than  those  for  a  R^  of  400.  Tliis  would  indicate 
the  presence  of  a  maximum  value  of  slip  velocity  somewhere  between  the 


1/3 


two  values. 

These  two  figures  eerve  to  illustrate  the  interrelationship 
between  holdup  (or  slip  velocity)  and  flow  pattern^  i.  e.  to  determine 
the  slip  velocity  occurring  in  the  vertical  flow  of  two  phase  mixtures 
one  must  first  determine  the  flow  pattern  present  in  the  two  pha.se  mixture. 
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SUPERFICIAL  WATER  VELOCITY,  V,  ,  ft.  per  sec. 
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TO 

CHAPTER  VII 

CONCLUSIONS 

(1)  Tests  have  heen  conducted  on  the  vertical  flow  of  air- 
water  mixtures  over  the  following  range  of  variables:  Tube  diameter 
0.63;  1*025;  1.50  and  2.5O  inches;  test  section  raid  point  pressure  3^ 
psia;  water  rates  from  .000149  to  O.O7OO  cubic  feet  per  second  and  gas 
rates  from  .00397  O.5IO  cubic  feet  per  second  at  average  test  section 
conditions. 

(2)  The  general  slope  of  the  unit  pressure  drop  curve  and  of 
the  holdup  curve  is  unchanged  with  varying  diameter.  However;  the 
difference  in  pressure  drop  between  the  maximum  and  two  minima  decreases 
with  increasing  diameter.  This  would  indicate  that  the  curve  would  be 
roughly  parabolic  in  shape  at  some  sufficiently  large  diameter. 

(3)  The  superficial  friction  factor;  f^  ;  has  been  successfully 
correlated  with  superficial  air  and  water  velocities  and  tube  diameter. 

(4)  Using  Figure  27  and  equ8>tion  2  and  4  one  may  calculate  the 
total  pressure  drop  for  the  vertical  flow  of  air-water  mix-cures  at  an 
average  pressure  of  36  psia  and  in  a  range  of  tube  diameters .between 
0.630  and  2.50  inches. 

(5)  The  flow  patterns  are  the  same  in  the  tube  diameters  studied 
with  o-nly  'bheir  transition  loci  and  the  mechanism  of  transition  from 

one  type  to  another  changing.  The  date  for  the  2.50  inch  tube  suggests 
that  the  froth  flow  pattern  may  disappear  entirely  at  some  larger 
diameter . 

(6)  The  holdup  ratio  for  the  vertical  flow  of  air -water  mixtures 
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may  be  computed  using  either  Figure  36  or  Figure  37  knowing  the  type 
of  flow  pattern.  The  latter  may  be  determined  using  Figures  33^  3^!- 
35*  luiowing  the  value  of  the  holdup  ratio,,  one  may  compute  the  actual 
slip  velocity  using  equation  7* 

(7)  If  more  data  were  available  at  very  low  air  rates  the 
manner  in  which  the  two  phase  friction  factor  becomes  coincident  with 
the  single  phase  friction  factor  could  be  studied. 

(8)  V/hen  the  effect  of  gas  phase  density  is  known,,  the  general 
method  applied  here  could  be  used  to  predict  the  pressure  drop  and  holdup 
for  the  vertical  flow  of  air-water  mixtures  under  any  conditions.  If  the 
variation  in  pressure  is  too  large  to  justify  the  use  of  an  average  value 
a  stepwise  procedure  could  be  employed. 
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APPEHDIX  A 

Table  AI  contains  additional  data  giving  the  lengths  and 
diameters  of  the  test  sections  and  the  volume  of  each  section  bet\7een 
plug  valves. 

The  basic  experimental  data  are  tabulated  in  Table  All.  The 
test  number  is  given  in  column  1  and  the  water  rate  in  cubic  feet  per 
second  in  column  2.  Colimin  3  reports  the  saturated  air  rate  in  cubic 
feet  per  second^  at  average  test  section  conditions.  Column  k  lists 
the  average  test  section  temperature.  The  discharge  air  water  ratio^ 
obtained  by  dividing  column  3  by  column  2j  is  given  in  column  5-  The 
Test  Section  air  water  ratio  is  reported  in  column  6.  Colimm  7  presents 
the  holdup  ratio^  being  the  ratio  of  column  5  "to  column  6.  The  pressure 
drop  mes.sured  in  feet  of  water  per  foot  of  tube  is  given  in  column  8 
and  the  flow  pattern  regime  reported  in  column  9* 
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TABLE  AI 

SUPPLEMENTARY  DATA 


Tube 

Diameter 

Test  Section 

Len{rbh 

Volume  Between 
Plug  Valves 

0.630  ins. 

52.88  ft. 

1940 

1.025  ins. 

52.88  ft. 

5290  c.c. 

1.50  ins. 

52.88  ft. 

10950 

2.50  ins. 

22.88  ft. 

28300 
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APPENDIX  B 


Tlie  results  calculated  from  Table  All  are  given  in  Table  BI 

of  Appendix  B.  Colunm  1  gives  the  test  number  and  the  unit  pressure 

drop  measured  in  feet  of  water  per  foot  of  tube  is  listed  in  column  2. 

The  hydrostatic  head  component^  1  Rm  and  the  irreversibility  com- 

1  +  Rv 

nonent^  _ 1 _ ,  (column  2  less  column  3)  both  measured 

1  -h  R,,  AXl 

in  feet  of  water  per  foot  of  tube^  are  reported  in  columns  3  and  4 
respectively.  The  superficial  water  velocity  in  feet  per  second  is 
given  in  coluim  5  while  the  superficial  gas  velocity^  also  in  feet  per 
second  is  listed  in  column  6.  Column  7  reports  the  valiie  of  the  super¬ 
ficial  friction  factor  and  the  product  DV^  is  given  in  column  8. 
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APPENDIX  C 

Table  Cl  contains  data  for  the  calibration  of  the  water 
rotameters.  The  rotameter  setting  is  given  in  coliimn  1.  The  mass  of 
water  collected  and  the  time  interval  is  shown  in  col-umns  2  and  3 
respectively.  The  calculated  water  rate  in  pounds  per  second  or 
cubic  feet  per  second  is  listed  in  column  4  and  column  5  reports 
the  value  of  W/H^  the  ratio  of  the  flow  rate  to  the  rotameter  setting. 
Plots  of  W/H  versus  rotameter  setting  are  given  in  Figures  AI  to  A4 
following  Table  1.  This  plot  was  chosen  rather  than  the  more  common 
plot  of  V/  versus  H  because  it  is  much  more  accurate. 

Table  CII  presents  data  for  the  calibration  of  the  orifice 
"Opiate  constants^  using  a  critical  flow  prover.  Column  1  shows  the 
air  orifice  plate  number  and/or  the  diameter  in  inches.  Column  2  reports 
the  pressure  upstream  of  the  air  orifice  plate  in  psia  and  the  manometer 
differential  in  centimeters  of  mercury  under  air  is  listed  in  column  3* 
(Some  of  the. latter  quantities  were  actually  read  as  centimetres  of 
water  under  air  but  have  been  converted  to  centimetres  of  mercury  for 
convenience).  Tlie  air  temperature- uiDstream  of  the  critical  flow  prover 
is  given  in  col-umn  4.  Column  5  presents  the  pressure  upstream  of  the 
critical  flow  prover  in  psia^  and  the  diameter  of  the  critical  flow 
orifice  is  listed  in  column  6.  The  value  of  the  critical  flow  constant 
for  the  orifice  given  in  column  6  is  listed  in  column  7* 

The  results  calcula^ted  from  the  data  in  Table  CII  are  given 
in  Table  CIII.  Column  1  is  the  same  for  both  tables.  Column  2  reports 
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the  flow  rate^  Q_,  in  cubic  feet  per  second.  The  value  of  Q  was  determined 


using  the  equation: 

(13): 

Q 

=  CP 

where  Q 

=  flow  rate^  1000  ft.'^/day  or  ft.'^/sec. 

at  l4,4  psia  and  60°F. 

C 

=  critical  flow  prover  coefficient 

corresponding  to  the  units  of  Q 

(The  numerical  values  for  C  were  obtained  from 

reference  (l3)* 

P 

=  pressure  upstream  of  the  critical  flow  prover,  psia 

S 

=  specific  gravity  of  the  gas  (l  for  air) 

T 

=  absolute  air  temperature  upstreajn  of  the 

critical  flow  prover, 

Z 

=  compressibility  factor 

The  numerical  value  of  the  term  Blh  calculated  from  columns 

T 

2.,  3  and  4  of  Table  II  is  listed  in  column  3.  Column  k  reports  the  value 
of  the  orifice  plate  constant  K.  It  was  calculated  using  the  equation 


Q 

K 

T 

vfnere  Q, 

3  ° 

=  air  flow  rate,  ft.  /sec.  at  lk»k  psia  and  60  F 

K 

=  air  orifice  plate  constant 

P 

=  pressure  upstream  of  air  orifice  plate,  psia 

=  me^nometer  differential^  inches  mercury  wonder  air 
T  =  upstream  air  temperature^ 


Figure  C5  following  Table  CIV  shows  a  plot  of  the  average 


value  of  K  for  each  orifice  plate  versus  the  orifice  diameter  in  inches. 
This  graph  serves  as  a  check  on  the  reliability  of  the  calibrations^ 
since  if  the  constants  are  accurate  they  should  all  line  on  a  straight 
line  with  a  slope  of  approximate!;^^  two. 

Table  CIV  lists  the  value  of  the  orifice  constant  E;  defined 
above ^  to  be  used  when  the  manometer  differential  (Z^)  is  in  centimeters 


of  mercury. 
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K  0  II  E  H  C  L  A  T  U  R  E 


C  = 

Cj)  = 

D  = 

4  = 

Z\F 

ZiF 


critical  flow  constant 
=  drag  coefficient 
^  diameter  of  pipe^  Ft. 

:  superficial  friction  factor  (tased  on  liquid  velocity) 
:  superficial  friction  factor  ("based  on  gas  velocity) 

=  VL 


gcD 


G  = 
Sc  ^ 

All 
K  : 
L  : 
Le 

hi 

P  = 
Alte 

Q  ^ 


1  2 
=  2fg  Vq 


gcD 


satiu'ated  air  rate^  To. /sec. 
jdimensional  conversion  faQtorj.  H  Aj'-  i't 


Ib.^  sec.^ 

=  holdup  ratio  (see  text) 

=  differential  manometer  reading;  inches  mercury  under  air 
:  air  orifice  constant 

-  waiter  rate;  lb .  /  sec . 

=  equivalent  leng-fch 

=  gas  liquid  mass  ratiO;  G/L 

=  gas  liquid  volume  ratiO;  GvG 

Lvl 

2  /  2 

-  pressure;  To. /ft.  absolute  or  lb/ in.  absolute 

loj.  2 

=  pressure  loss  due  to  drag  (  Ih/ft.  ) 

=  rate  of  gas  floW;  ft. •^/ sec.  or  1000  ft. day 
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S  =  specific  gravity  of  the  gas 

T  =  Temperatiire  ^  ^ 

=  gas  phase  velocity,  based  on  average  flow  conditions  and 
the  empty  cross  section  of  the  tube 
Yj^  =  liquid  velocity  based  on  the  empty  cross  section  of  the 
tube 

Vq  =  average  lineal velocity  of  gas  based  on  fraction  of  tube 
filled  with  ,gas 

vi  =  average  lineal  velocity  of  liquid  based  on  fraction  of 
tube  filled  with  liquid 
X  =  tube  heiglit,  ft.  ‘ 

Z  =  compressibility  factor 

3 

pQ  =  density  of  gas  at  flowing  conditions,  lb./ ft, 

=  density  of  liquid,  lb/ ft. 

[1^  =  viscosity  of  liquid,  lb. m/ft. sec. 

[1^  =  viscosity  of  gas  at  flowing  conditions,  Ib.m/ft.sec. 

3 

Vj^  =  specific  volimie  of  liquid,  ft.  /lb. 

3  / 

Vq  =  specific  volume  of  gas  at  flowing  conditions,  ft.  /lb. 

^  =  sLU’face  tension 


5 


tube  roughness 
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